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Abstract 


x  In  this  paper,  we  consider  a  multicommodity  flow  problem  where  for  each  pair  of  vertices, 

( u,v ),  we  are  required  to  send// half-units  of  commodity  (u,v)  from  u  to  v  and  /  half-units  of 
commodity  (v,u)  from  v  to  u  without  violating  capacity  constraints.  Our  main  result  is  an 
algorithm  for  performing  theltask  provided  that  the  capacity  of  each  cut  exceeds  the 
demand  across  the  cut  by  a  6(log  n)  factor.  The  condition  on  cuts  is  required  in  the  worst 
case,  and  is  trivially  within  a  e>(log  n}  factor  of  optimal  for  any  flow  problem. 


The  result  is  of  interest  because  it  can  be  used  to  construct  the  first  poly-log  times  optimal 
approximation  algorithms  for  a  wide  variety  of  problems,  including  minimum  quotient 
separators,  1/3-2/3  separators,  bifurcators,  crossing  number  and  VLSI  layout  area.  The 
result  can  also  be  used  to  efficiently  route  packets  in  arbitrary  distributed  networks.  For 
example,  we  can  prove  that  any  n-node  bounded  degree  graph,  G,  with  minimum  edge 
expansion  fa  can  be  configured  off-line  to  simulate  any  w-node  bounded  degree  graph  H  in 
0(log  n/a)  steps  using  constant  size  queues.  By  letting  H  be  a  universal  network,  we  can 
then  use  G  to  simulate  a  PRAM  on-line  with  delay  0(log2  nj  a)  steps,  resolving  a 
conjecture  in  the  theory  of  distributed  computation.  These  results  compare  well  with  the 
universal  lower  bound  of  n(l/a)  and  the  existential  lower  bound  of  o(log  n/a). 
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Abstract 

In  this  paper,  we  consider  a  multicommodity 
flow  problem  where  for  each  pair  of  vertices, 
(u,  v),  we  are  required  to  send  /  half-units  of 
commodity  (ti,i>)  from  u  to  c  and  /  half-units 
of  commodity  (v.u)  from  v  to  u  without  vio¬ 
lating  capacity  constraints.  Our  main  result  is 
an  algorithm  for  performing  the  task  provided 
that  the  capacity  of  each  cut  exceeds  the  de¬ 
mand  across  the  cut  by  a  0(logn)  factor.  The 
condition  on  cuts  is  required  in  the  worst  case, 
and  is  trivially  within  a  0(logn)  factor  of  op¬ 
timal  for  any  flow  problem. 

The  result  is  of  interest  because  it  can  be 
used  to  construct  the  first  poly-log  times  opti¬ 
mal  approximation  algorithms  for  a  wide  vari¬ 
ety  of  problems,  including  minimum  quotient 
separators,  3-5  separators,  bifurcators,  cross¬ 
ing  number  and  VLSI  layout  area.  The  result 
can  also  be  used  to  efficiently  route  packets  in 
arbitrary  distributed  networks.  For  example, 
we  can  prove  that  any  n-node  bounded  degree 
graph,  G ,  with  minimum  edge  expansion  a  can 
be  configured  off-line  to  simulate  any  n-node 
bounded  degree  graph  H  in  0(^2)  steps  us¬ 
ing  constant  size  queues.  By  letting  H  be  a 
universal  network,  we  can  then  use  G  to  simu¬ 
late  a  PRAM  on-line  with  delay  0(^-a)  steps, 
resolving  a  conjecture  in  the  thoery  of  dis¬ 
tributed  computation.  These  results  compare 
well  with  the  universal  lower  bound  of  fi(£) 
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and  the  existential  lower  bound  of 

1  Introduction 

1.1  Multicommodity  Flow  Prob¬ 
lems 

A  multicommodity  flow  problem  (MFP)  con¬ 
sists  of  a  graph  G  =  (V',£),  a  set  of  com¬ 
modities,  T,  a  capacity  function  on  the  edges 
C  :  E  -*  R+ ,  a  supply  function  for  vertices 
5  :  V  x  C  -*•  R+,  and  a  demand  function  for 
vertices  D  :  V  x  T  -*  R+ .  The  object  is  to  con¬ 
struct  flows  for  the  commodities  that  satisfy 
the  demand  for  each  commodity  at  each  vertex 
without  violating  the  constraints  imposed  by 
the  supply  function  and  the  capacity  function. 

It  is  well  known  that  any  MFP  can  be  solved 
in  polynomial  time  using  linear  programming. 

In  this  paper,  we  consider  a  special  kind  of 
MFP  that  we  call  a  uniform  multicommodity 
flow  problem  (UMFP).  In  a  UMFP,  we  have 
a  separate  commodity  for  each  ordered  pair  of 
vertices  (i.e.,  r  =  V  x  V),  and  we  insist  that  / 
half- units  of  commodity  (u,t>)  flow  from  u  to  v 
and  that  /  half-units  of  commodity  (v,  u)  flow 

This  research  *u  supported  in  put  by  the  De¬ 
fense  Advanced  Research  Projects  Agency  under  Con¬ 
tract  N 0001 4-80-00622  the  Office  of  Naval  Research 
under  Contract  N00014-86-K-0593,  the  Air  Force  un¬ 
der  Contract  OSR-86-0076,  and  the  Army  under  Con¬ 
tract  DAAL-03-86-K-0171.  Tom  Leighton  is  also  sup-  / 
ported  in  put  by  an  NSF  Presidential  Young  Investiga-  1 
tor  Awud  with  matching  funds  from  IBM  and  AT&T.  .. 


from  v  to  u  for  each  u.r  €  V.  The  underly¬ 
ing  graph  and  capacities  can  be  arbitrary,  al¬ 
though  we  will  only  need  the  special  case  when 
the  edge  capacities  are  all  equal  for  the  appli¬ 
cation.  The  object  of  a  UMFP  is  to  maximize 
/  subject  to  the  constraints  imposed  by  the 
edge  capacities.  This  can  be  done  in  polyno¬ 
mial  time  using  linear  programming,  and  can 
be  approximately  done  using  faster  algorithms 
by  Shahrokhi  and  Matula  [19]. 

As  in  the  case  for  the  single  commodity 
flows,  an  MFP  cannot  be  solved  if  there  is  a  cut 
of  the  graph  for  which  the  demand  across  the 
cut  exceeds  the  capacity  of  the  cut.  In  the  case 
of  a  UMFP,  the  demand  across  a  cut  (U,~U)  is 
simply  /  |r|  |F|  where  U  C  V.  The  capacity 
of  the  cut  is  simply  the  sum  of  the  capacities 
of  the  edges  crossing  the  cut.  In  the  case  when 
all  the  capacities  are  equal,  the  capacity  of  the 
cut  is  proportional  to  the  number  of  edges  in 
the  cut.  This  motivates  the  definition  of  the 
minimum  cut  for  a  UMFP  to  be 


Note  that  for  any  feasible  solution  to  a  UMFP, 
we  must  have  f  <  S. 

The  value  of  S  defined  in  equation  1  also  has 
interpretations  in  graph  theory.  If  C(e)  =  1 
for  all  e  6  E,  then  S  is  the  value  of  the  spars¬ 
est  cut  of  the  graph.  In  the  case  of  arbitrary 
edge  capacities,  then  S  is  the  weighted  spars¬ 
est  cut  value  for  the  graph  with  edge  weights 
given  by  C.  Even  in  the  simple  case  where  all 
edge  weights  (capacities)  are  1,  computing  the 
sparsest  cut  of  a  graph  is  N  P- hard  [14]. 

A  closely  related  quantity  to  the  sparsest  cut 
of  a  graph  is  the  minimum  edge  expansion  or 
flux  of  a  graph,  defined  by 


connected  to  the  rest  of  the  graph  with  edges 
of  total  weight  at  least  a  \U\.  For  example,  un¬ 
weighted  (i.e.,  C(e)  =  1  for  all  e  e  E)  expander 
graphs  have  constant  flux.  A  cut  that  achieves 
the  flux  is  called  the  minimum  quotient  sepa¬ 
rator  and  is  related  by  a  constant  factor  to  the 
sparsest  cut  since 

| S<a<nS  (3) 

for  any  n-vertex  graph  since  the  number  of  ver¬ 
tices  in  the  big  half  of  a  cut  must  lie  between 
§  and  n.  Computing  the  minimum  quotient 
separator  is  also  ATP-complete.  (The  proof  is 
identical  to  that  of  graph  bisection  [7].) 

Unlike  the  case  for  single  commodity  flow 
problems,  there  is  no  max-flow,  min-cut  the¬ 
orem  for  the  UMFP.  In  particular,  the  condi¬ 
tion  that  S  >  f  is  necessary  but  not  always 
sufficient  for  a  UMFP  to  have  a  flow  of  size 
/.  As  a  particularly  nasty  example,  consider 
a  UMFP  on  a  bounded-degree  expander  graph 
with  identical  edge  capacities  C(e )  =  C0.  Since 
most  pairs  of  vertices  are  separated  by  distance 
fi(logn)  in  such  a  graph,  the  sum  of  the  edge 
capacities  must  be  n(/n2logn;  in  order  for 
there  to  be  a  flow  of  size  /  in  the  UMFP.  This 
means  that  Co  =  fi(/nlogn)  which  is  equiva¬ 
lent  to  requiring  that  S  =  fi(/logn).  Hence, 
in  this  case,  the  capacity  across  each  cut  must 
exceed  the  demand  across  the  cut  by  a  0(log  n) 
factor  in  order  for  there  to  be  a  solution  to  the 
UMFP. 

The  preceding  example  is  as  bad  as  things 
ca.  «;et  however.  This  fact  is  summarized  in 
the  ang  theorem. 

Theorem  1  There  is  a  feasible  flow  with  f  = 
^(iofn)  for  any  ^-vertex  UMFP. 

The  preceding  result  provides  an  ap¬ 
proximate  max-flow  min-cut  relationship  for 
UMFPs  that  is  tight  in  the  worst  case.  In  par¬ 
ticular,  we  have  that 


In  other  words,  a  graph  has  flux  a  if  every 
subset  U  with  at  most  half  of  the  vertices  is 


(4) 


* 


for  any  UMFP.  In  other  words,  the  maximum 
flow  is  always  within  a  0(logn)  factor  of  the 
minimum  cut.  Moreover,  there  are  examples 
of  both  /  =  0(5 )  and  /  =  0(5/  log  n),  so  the 
approximation  is  existentially  tight. 

The  proof  of  Theorem  1  is  the  first  of  two 
main  results  in  the  paper  and  is  included  in 
Section  2.  In  the  remainder  of  this  section, 
we  describe  the  method  of  the  proof,  the  sec¬ 
ond  main  result,  and  the  applications  of  the 
results  to  a  wide  variety  of  problems.  Section 
3  contains  the  proof  of  Theorem  2.  We  con¬ 
clude  with  some  remarks,  acknowledgements, 
and  references  in  Sections  4  -  6. 

1.2  Proof  Method 


proving  Theorem  1  is  showing  that  paths  with 
small  total  weight  exist  no  matter  how  the  ad¬ 
versary  places  the  weights.  It  is  also  necessary 
for  some  of  the  applications  to  find  the  shortest 
paths  in  a  way  that  reveals  a  cut  whose  size  is 
inversely  proportional  to  the  total  path  length. 
Hence  we  must  use  a  nonstandard  short  paths 
type  of  algorithm. 

We  suspect  that  our  method  for  proving  the 
existence  of  paths  with  short  total  length  may 
be  directly  useful  in  other  contexts  such  as 
routing  paths  around  faults  or  costly  (e.g.,  con¬ 
gested)  components  in  a  distributed  network. 
In  this  paper,  however,  we  use  this  technique 
only  indirectly  -  in  the  context  of  the  dual 
problem  of  finding  sparsest  cuts. 


To  prove  Theorem  1,  we  work  with  the  dual 
of  a  UMFP.  The  dual  UMFP  problem  consists 
of  a  weighted  graph  where  the  weight  of  an 
edge  (denoted  by  d(e))  is  treated  as  a  distance 
and  is  selected  by  the  adversary  subject  to  the 
condition  that  £*e£:  d(e)C{e)  =  1.  The  object 
is  to  route  the  commodities  along  paths  in  a 
way  that  minimizes  the  total  distance  of  all 
paths  used.  In  general,  a  UMFP  has  a  feasible 
flow  /  if  and  only  if  for  every  legal  distance 
function,  the  sum  of  all  path  distances  is  at 
most  j.  (This  fact  can  be  proved  using  the 
duality  theory  of  linear  programming  [9].)  In 
other  words. 


/  = 


max**  {minimum  total  path  length} 


.(5) 


Given  an  n-node  graph  with  capacities  C 
and  weighted  flux  a,  and  any  distance  func¬ 
tion  d  such  that  J y€E  d(e)C(e)  =  1,  we  show 
in  Section  2  how  to  route  paths  between  all 
pairs  of  vertices  so  that  the  total  path  weight 
is  at  most  0(n-^p )  =  This  will  be 

sufficient  to  prove  Theorem  1  since  we  can  then 
conclude  that  there  is  a  feasible  flow  of  size 

It  is  worth  noting  that  it  is,  of  course,  well- 
known  how  to  compute  shortest  paths  in  a  net- 
v,*ork  in  polynomial  time.  The  tricky  part  of 


1.3  Approximating  Flux 

From  Equations  3  and  4,  we  can  quickly  deduce 
that 

Y  <  a  <  0(/nlogn)  (6) 

for  any  weighted  gTaph.  Hence,  we  can  approx 
imate  the  flux  of  any  graph  up  to  a  0(logn)- 
factor  simply  by  finding  the  optimal  value  of  / 
for  the  associated  UMFP.  Unfortunately,  this 
approach  is  not  sufficient  to  actually  locate 
an  approximate  minimum  quotient  separator, 
however. 

In  order  to  find  am  actual  cut,  we  first  solve 
the  corresponding  UMFP  and  its  dual  in  order 
to  construct  the  worst  possible  distance  func¬ 
tion  d  :  E  — »  R+ .  We  then  use  the  path  routing 
algorithm  described  in  Section  2  to  find  paths 
between  all  pairs  of  vertices  with  small  total 
path  length.  As  the  algorithm  progresses,  we 
keep  track  of  the  minimum  edge  expansion  cut 
found  up  to  that  point.  At  the  end,  we  out¬ 
put  the  best  such  cut  ever  encountered.  The 
algorithm  has  the  nice  property  that  the  total 
path  length  is  at  most  where  a is 

the  sparsest  cut  observed  during  the  running 
of  the  algorithm.  Hence  by  Equations  5  and  6 
we  have 


Qobs  <  O(otlogn), 

and  thus  our  minimum  observed  sparsest  cut 
is  within  a  0(logn)  factor  of  optimal. 

1.4  Approximating  Separators,  Bi- 
furcators,  Crossing  Numbers 
and  VLSI  Layout  Area 

The  algorithm  for  finding  an  approximate  min¬ 
imum  quotient  cut  can  be  used  to  find  polylog 
approximation  algorithms  for  a  variety  of  N P- 
hard  combinatorial  problems.  For  example, 
bv  recursively  applying  the  approximate  mini¬ 
mum  quotient  algorithm  as  in  [18],  we  can  ob¬ 
tain  a  0(log2n)  times  optimal  approximation 
algorithm  for  finding  minimum  graph 

separators  for  any  k  >  3.  Alternatively,  we 
can  obtian  obtain  a  O(log  n)  times  optimal  ap¬ 
proximation  algorithm  for  finding  approximate 
£,  separators.  (The  difference  in  the  lat¬ 
ter  result  is  that  we  pay  the  price  for  a  £ , 
separator  but  only  get  a  separator  as 

output  where  a  >  1  is  a  small  constant.) 

These  algorithms  can  then  be  combined  with 
the  reductions  of  Bhatt  and  Leighton  [1]  to  ob¬ 
tain: 

1.  a  0(log4n)  times  optimal  approximation 
algorithm  for  c  +  n  where  c  is  the  cross¬ 
ing  number  of  an  n-node  bounded  degree 
graph, 

2.  a  0{ log5/2  n)  times  optimal  approxima¬ 
tion  algorithm  for  the  >/2-bifurcators  of 
an  n-node  bounded  degree  graph,  and 

3.  a  O(log6n)  times  optimal  approximation 
algorithm  for  the  layout  area  of  an  n-node 
bounded  degree  graph. 

At  first  glance,  the  crossing  number  approx¬ 
imation  algorithm  seems  fairly  weak  since  we 
are  only  approximating  c  +  n  and  only  for 
bounded  degree  graphs.  However,  c  >  w(n) 
for  any  graph  with  bisection  width  *>(y/n)  as 


well  as  any  graph  with  4n  or  more  edges  [11]. 
These  conditions  are  satisfied  by  most  graphs. 
Moreover,  the  result  can  be  partially  extended 
to  graphs  with  arbitrary  degree.  For  example, 
the  same  result  holds  for  regular  (or  near  reg¬ 
ular)  graphs  of  arbitrary  degree. 

We  should  note  that  some  minor  effort  is  re¬ 
quired  to  adapt  the  arguments  of  Bhatt  and 
Leighton  [1]  to  obtain  these  results.  Full  de¬ 
scriptions  of  the  algorithms  and  their  analysis 
will  be  included  in  the  full  paper. 

1.5  Application  to  Packet  Routing 
in  Distributed  Networks 

Communication  in  a  distributed  network  is  an 
important  area  of  research  for  which  relatively 
little  is  known.  For  example,  given  an  arbitrary 
n-processor  network,  we  might  wish  to  know 
how  well  it  can  simulate  a  well-studied  network 
or  other  parallel  machine  such  as  the  butterfly 
or  PRAM.  Some  progress  has  been  made  on 
this  problem  in  the  special  case  that  the  graph 
is  an  expander.  For  example.  Peleg  and  Upfal 
have  shown  how  to  solve  any  n-packet  routing 
problem  on  any  n-node  expander  in  O(logn) 
expected  steps  using  queues  of  size  O(logn)  at 
each  node  [15].  For  general  graphs,  relatively 
little  is  known  except  that  it  has  been  observed 
by  many  researchers  that  flux  and  diameter 
are  important  parameters  that  influence  per¬ 
formance  (e.g.,  they  are  both  lower  bounds  on 
the  time  needed  to  route  a  random  permuta¬ 
tion). 

In  this  paper,  we  provide  the  first  step 
towards  a  general  solution  to  the  problem 
by  showing  how  any  n-node  bounded-degree 
graph,  G,  with  flux  a  can  simulate  any  other 
n-node  bounded-degree  graph,  ff,  with  delay 
0(^5^)  constant  size  queues.  The  result 
is  optimal  in  the  sense  that  there  are  simula¬ 
tion  problems  that  require  this  much  time,  and 
is  robust  in  the  sense  that  the  simulation  can 
take  place  for  any  1-1  embedding  of  the  nodes 
of  H  onto  nodes  of  G.  The  main  drawback  of 


the  simulation  result  is  that  it  requires  off-line 
computation.  However,  once  the  off-line  em¬ 
bedding  is  performed  for  one  H,  simulation  of 
another  H'  can  be  performed  on-line  by  using 
H  to  simulate  H'.  For  example,  by  embed¬ 
ding  a  butterfly  or  other  uni  versal  network  into 
G  off-line,  we  can  then  use  G  to  simulate  any 
CRCW  PRAM  algorithm  in  an  on-line  fashion 
with  delay  0(  lo^-~ )  and  constant  size  queues. 

Conceivably,  such  a  result  could  have  a  sub¬ 
stantive  impact  on  the  theory  of  distributed 
computation,  where  the  complexity  of  algo¬ 
rithms  for  operations  such  as  sorting  is  typi¬ 
cally  measured  in  terms  of  numbers  of  nodes 
and  edges  and  (sometimes)  diameter.  While 
such  algorithms  are  of  interest  lor  some  graphs, 
they  can  be  far  from  optimal  for  many  others. 
By  using  our  flux-based  approach,  however, 
it  is  possible  to  devise  a  single  on-line  algo¬ 
rithm  to  sort  in  any  network  in  O(log2  n)  steps 
times  the  optimal  for  that  graph.  Of  course, 
we  will  still  need  off-line  computation  to  set 
up  routing  tables  for  the  graph  (in  this  case, 
we  embed  the  AKS-based  sorting  network  de¬ 
scribed  by  Leighton  in  [13]),  but  if  the  sorting 
is  to  be  done  many  times,  the  algorithm  will 
be  asymptotically  much  better  than  currently 
known  techniques  for  general  graphs.  (Alter¬ 
natively,  we  could  use  0(log3  n)  times  optimal 
steps  and  obtain  a  smaller  constant  factor  by 
off-line  embedding  a  shuffle-exchange  graph.) 

In  order  to  prove  that  any  n-node  bounded- 
degree  graph  G  with  flux  a  can  simulate  any 
other  n-node  bounded-degree  graph  H  with 
delay  0(^1),  we  first  prove  that  H  can  be 
embedded  in  G  with  congestion  and  dilation 
<9(^p).  This  is  Theorem  2  and  is  the  second 
of  our  two  main  results.  (The  congestion  of  an 
embedding  is  the  maximum  number  of  edges  of 
H  that  traverse  any  edge  of  G.  The  dilation  of 
an  embedding  is  the  maximum  distance  sepa¬ 
rating  adjacent  nodes  of  H  when  embedded  in 
G .)  The  proof  is  then  completed  by  applying 
the  recent  result  of  Leighton,  Maggs  and  Rao 
[12]  that  any  packet  routing  problem  can  be 


solved  in  0( congestion  +  dilation)  steps  using 
constant  size  queues. 

Unfortunately,  the  LMR  routing  algorithm 
is  nonconstructive.  We  can  circumvent  this 
drawback  by  introducing  randomness  and 
queues,  however.  In  particular,  we  can  use 
Theorem  2  to  embed  Kn  in  G  with  congestion 
0(5i2£IL)  and  dilation  O(^).  This  can  be 
performed  in  polynomial  time  off-line.  We  can 
then  route  a  random  permuation  in  0(  - ) 

expected  steps  using  queues  of  size  )  by 

introducing  a  random  delay  selected  uniformly 
from  [l,^p]  for  each  packet  and  multiplexing 
the  packets  as  they  arrive  at  each  edge.  The 
details  are  not  difficult  to  work  out  and  will  be 
included  in  the  final  paper. 

One  final  comment  is  relevant  here.  The 
observant  reader  will  notice  that  we  are  rout¬ 
ing  arbitrary  paths  in  G  within  time 
without  regard  for  the  diameter  of  G.  At  first 
glance,  this  would  not  seem  to  be  possible  if 
the  diameter  were  u/(i2&H).  Such  a  scenario  is 
not  possible,  however,  since  as  a  consequence 
of  our  proof  we  will  find  that  the  diamete-  of 
any  bounded-degree  graph,  is  always  0(^£p). 

1.6  Previous  Work 

Aside  from  the  expander- related  results  cited 
in  Section  1.5,  there  axe  relatively  few  previous 
results  (that  we  are  aware  of)  on  the  problems 
addressed  in  this  paper.  Although  there  are 
many  max-flow  min-cut  results  for  multicom¬ 
modity  flow  problems  [14]  [6],  they  tend  to  be 
exact  results  (or  sometimes  “plus  one”  results) 
for  specialized  problems,  and  we  know  of  noth¬ 
ing  resembling  the  approximate  result  in  this 
paper.  Duals  to  max-flow  problems  that  in¬ 
volve  min  cuts  have  been  constructed  (see  [19]), 
but  they  involve  distance  functions  and  are  not 
natural  min-cut  problems  in  the  normal  sense. 

Although  there  are  numerous  heuristic  algo¬ 
rithms  for  separator  related  problems  [5][10], 
as  well  as  numerous  average-case  algorithms  [4] 
[2]  [3],  there  a  very  few  previous  results  in  ap- 


proximation  algorithms.  For  example,  the  best 
previously  known  algorithms  for  separators  are 
due  to  Plaisted  [16],  who  has  a  randomized 
0(v/n)  times  optimal  approximation  algorithm 
and  a  randomized  O(logn)  times  the  square  of 
optimal  approximation  algorithm.  There  has 
been  substantial  recent  progress  on  the  related 
problem  of  bandwidth  minimization,  however. 
In  particular,  Gonzalez  and  Kurki-Gowdara 
have  used  bundle-capacity  flow  techniques  to 
construct  a  2  times  optimal  approximation  al¬ 
gorithm  for  this  problem  as  well  as  for  a  vari¬ 
ety  of  VLSI  layout  problems  in  [8].  Whether 
or  not  these  techniques  can  be  combined  with 
those  presented  in  this  paper  to  obtain  better 
approximation  algorithms  for  separator  related 
problems  is  a  potentially  interesting  area  of  re¬ 
search. 


2  Proof  of  Theorem  1 

Following  the  analysis  of  Section  1.2,  we  need 
only  show  that  for  any  n- node  graph  G  - 
(V,  E)  with  weights  C  :  E  —  R+  and  dis¬ 
tance  function  d  :  E  — *  R+  such  that 
£e€£  d(e)C{e)  =  1,  it  is  possible  to  route 
paths  between  every  pair  of  vertices  with  total 
distance  0(nl°gn )  where  a  is  the  weighted  flux 
of  G.  This  equivalent  to  showing  that  there  is 
some  vertex  u  of  G  for  which 


The  algorithm  for  finding  such  a  node  u  con¬ 
sists  of  three  Stages.  The  first  Stage  is  a  pre¬ 
processing  stage.  We  start  by  placing  an  in¬ 
tegral  number  of  “distance  tokens"  on  each 
edge.  In  particular,  we  place  max  jl,  j 
tokens  on  edge  e  for  every  e  €  E  where  s  = 
t-* — V,.  Each  token  corresponds  to  distance 

s  and  serves  to  discretize  the  distance  func¬ 
tion.  By  using  the  tokens,  we  can  break  up 
long  edges  into  discrete  pieces,  each  with  the 
same  amount,  s,  of  distance.  Note  that  the  to¬ 
tal  distance  of  each  edge  is  not  decreased  by 
this  process. 

We  next  assign  a  weight  to  each  token.  In 
particular,  tokens  on  edge  e  are  assigned  weight 
C(e).  Note  that  the  total  weight  of  all  the  to¬ 
kens  is  at  most 


e€E 


since  £e6£C(e)d(e)  =  1  and  ]Pe€£  C(e)  =  i. 

We  now  proceed  with  the  second  Stage  of  the 
algorithm.  The  purpose  of  Stage  2  is  to  grow  a 
shallow  depth  region  that  contains  a  majority 
of  the  nodes.  We  do  this  by  insisting  that  every 
time  we  grow  the  region  by  a  token,  we  make 
the  weight  of  the  region  heavier  by  an  amount 
that  is  proportional  to  the  expansion  ratio.  In 
particular,  we  execute  the  following  steps: 


where  d(u,v )  denotes  the  distance  along  the 
selected  path  from  u  to  v.  For  later  reference 
we  state  this  result  as  Lemma  1. 

Lemma  1  There  is  some  vertex  u  of  G  for 
which 

ffd(u,v)  =  0(^), 

where  d(u,v)  denotes  the  distance  along  the  se¬ 
lected  path  from  u  to  v. 


1.  Select  a  vertex  v  €  G,  to  initiate  the  re¬ 
gion. 

2.  Extend  the  region  by  one  token. 

3.  Let  i  denote  the  depth  (i.e.,  radius)  of  this 
region  as  measured  in  tokens.  (This  is  the 
number  of  times  we  have  executed  Step 
2.)  Let  Wi_j  denote  the  weight  of  tokens 
in  the  first  i  -  1  levels  of  the  region.  If 
Wi  >  (1  +  e)Wj_i  where  e  =  then  go 
to  Step  2. 

4.  If  the  region  contains  ^  or  more  nodes, 
halt  and  output  the  region. 


5.  Otherwise,  remove  the  nodes  in  the  region 
from  the  gtaph  and  go  to  Step  1. 


(anything  larger  could  never  be  used  in  a 
UMFP)  and  thus 


We  first  prove  that  the  algorithm  eventually 
halts  and  outputs  a  region  containing  at  least 
2  nodes.  For  the  purposes  of  contradiction,  as¬ 
sume  otherwise  and  consider  the  union  of  the 
removed  regions  at  the  first  instant  that  it  con¬ 
tains  at  least  7  nodes.  Since,  by  assumption, 
no  single  removed  piece  contains  ^  removed 
nodes,  the  union  of  the  removed  regions  has 
l  nodes  where  7  <  /  <  at  this  point.  By 
the  definition  of  a,  this  means  that  at  least  2p 
token  weight  leaves  the  union  of  the  regions. 
From  the  condition  at  Step  3,  we  know  that 
the  ratio  of  the  internal  token  weight  to  the 
boundary  token  weight  of  each  removed  piece 
is  greater  than  K  Hence  the  same  is  true  of  the 
union  of  the  removed  regions,  and  we  conclude 
that  the  union  of  the  removed  regions  contains 
at  least  77  total  token  weight.  This  is  impos¬ 
sible,  however,  since  77  =  JL  >  1  which  is  the 
total  token  weight.  Hence  the  algorithm  halts 
and  successfully  outputs  a  region  with  at  least 
^  nodes. 

We  next  show  that  the  depth  of  the  output 
region  (as  measured  in  numbers  of  tokens)  is 
at  most  0(^®-2-).  To  see  why,  we  need  only  ob¬ 
serve  that  W,  >  (1  +  €)W,_i  for  each  level  of 
the  output  region,  and  hence  that  a  region  of 
depth  d  must  have  total  token  weight  exceeding 
(1  +  «)<i“1cmm,  where  cmi„  is  the  minimum  to¬ 
tal  capacity  of  the  edges  incident  to  any  node. 
Since  there  is  only  7  =  2  C(«) total  token 

weight  we  can  conclude  that 

e€E 

(l  +  €)<‘-1cmin  <2£C(e), 
and  thus  that 


d  =  0( 


Jog 


T'*EC(e) 


cmn 


log(  1  +  () 


*€E 


52  ^  ”3c" 


SO 


log£— C(e)  =  0(logn). 


In  addition, 


nsa 


€  = 


Cmin 


<  -  1 

9£e€i  C(ej  “  9  < 


na 


since  by  averaging  over  all  nodes,  there  must  be 
some  single  node  for  which  the  weighted  min- 

2  C(e) 

imum  edge  expansion  is  at  most  ++— — — . 
Thus  log(l  +  e)  =  0(c),  and  d  =  0(J2£2)  as 
claimed. 

At  this  point,  we  have  constructed  a  region 
which  contains  at  least  ^  nodes,  and  for  which 
the  sum  of  the  distances  from  the  root  to  all 
nodes  in  the  region  is  at  mo6t  0(  —  -*  -)  = 
It  remains  to  extend  the  region  to  the 
remaining  nodes  with  short  total  path  length. 
This  takes  place  during  Stage  3. 

Stage  3  consists  of  simply  extending  the  re¬ 
gion  constructed  in  Stage  2  in  a  breadth  first 
manner  (one  level  of  tokens  at  a  time)  until  all 
nodes  are  reached.  Although  individual  paths 
constructed  during  Stage  3  can  be  much  longer 
than  those  constructed  in  Stage  2  (and  this  can 
be  required  in  certain  examples),  their  total 
length  will  still  be  small. 

To  prove  that  the  total  path  length  is  small 
for  the  full  region,  we  let  nj  denote  the  number 
of  nodes  still  not  in  the  region  after  it  is  grown 
by  t  tokens  during  Stage  3.  For  example,  no  is 
the  number  of  nodes  not  in  the  region  at  the 
end  of  Stage  2.  Since  precisely,  n,  paths  will 
traverse  the  tth  level  of  Stage  3,  the  total  path 
length  contributed  by  nodes  in  Stage  3  is  at 
most 


Without  loss  of  generality  we  can  assume 
that  the  capacity  of  any  edge  is  at  most  ncm,„ 


o(!ai)+.f 

i=0 


The  first  term  comes  from  the  portion  of  the 
paths  in  the  Stage  2  portion  of  the  region,  and 
the  second  term  comes  from  the  portion  of  the 
paths  routed  in  Stage  3. 

The  sum  n«  cannot  he  too  large,  how¬ 
ever,  since  if  after  level  t,  there  are  still  n,  un¬ 
seen  nodes,  the  weight  of  all  the  tokens  used  at 
level  i  must  be  at  least  an,.  Since  the  total  to¬ 
ken  weight  is  at  most  |,  we  can  conclude  that 


n,  <  and  thus  that 


i=0 


Hence  the  sum  of  all  path  lengths  is  at  most 
0(^p).  35  claimed. 

This  completes  the  proof  of  Lemma  1  and 
hence  Theorem  1.  In  order  to  construct  the  al¬ 
gorithm  for  finding  approximate  sparsest  cuts 
as  described  in  Section  1.3,  however,  we  still 
need  to  check  that  the  total  weight  of  the  paths 


in  the  region  is  at  most  where  aotj  is 


the  sparsest  cut  observed  during  the  construc¬ 
tion  of  the  region.  This  is  easily  done,  since  we 
simply  replace  a  with  a in  all  of  the  analy¬ 
sis  for  the  performance  of  the  algorithm.  Note 
the  algorithm  is  run  the  same  way  (c  doesn’t 
change)  but  the  analysis  can  be  improved  by 
using  a0b,  instead  of  a. 


3  Proof  of  Theorem  2 


Our  second  main  result  can  be  formally  stated 
as  follows. 


Theorem  2  Given  any  n-node  bounded  degree 
graph  H ,  and  any  1-1  embedding  of  the  nodes 
of  H  onto  the  nodes  of  G,  the  edges  of  H  can 
be  routed  as  paths  in  G  with  congestion  and 
dilation  0(12^). 


By  Lemma  1  and  standard  duality  argu¬ 
ments  we  can  deduce  that  the  commodities 
induced  by  the  edges  of  H  can  be  satisfied 
in  G  with  congestion  OC^1).  (Here  conges¬ 
tion  corresponds  to  capacity  in  the  UMFP.)  We 


can  then  make  use  of  Raghavan's  [17]  method 
of  converting  fractional  flows  into  single  paths 
to  show  that  H  can  be  embedded  in  G  with 
congestion.  These  peripheral  details 
will  appear  in  the  full  paper. 

In  order  to  produce  paths  with  simultaneous 
congestion  and  dilation  OC^2-),  however,  we 
have  to  extend  the  proof  of  Lemma  1  so  that  all 
paths  used  have  at  most  Ol^f2)  edges  in  ad¬ 
dition  to  total  (token)  distance  OC^2).  The 
result  will  then  follow  by  restricting  the  lin¬ 
ear  program  so  that  only  paths  with  distance 
are  allowed. 

To  construct  a  region  with  small  edge  depth 
as  well  as  small  total  path  distance  (measured 
as  in  Section  2  with  the  distance  function),  we 
need  to  modify  the  region  construction  algo¬ 
rithm  of  Lemma  1.  In  particular,  we  need  only 
modify  stage  3  since  the  edge  depth  of  the  re¬ 
gion  produced  by  Stage  2  is  bounded  by  the 
token  depth,  which  was  shown  to  be  at  most 
O(isp)  =  O(^).  This  is  at  most  0(!sp) 
since  s  =  ©(£)  for  an  n- node  bounded-degree 
graph  with  edge  capacities  all  of  which  are  one. 

To  modify  Stage  3,  we  examine  each  level  to 
see  whether  or  not  the  number  of  nodes  on  the 
level  exceeds  the  fraction  —  of  the  number 

*®mat 

of  tokens,  where  dmax  is  the  maximum  degree 
of  the  graph  G.  If  not,  then  we  insert  addi¬ 
tional  tokens  in  front  of  all  the  nodes  that  were 
on  the  level  in  order  to  “push”  them  to  the  next 
level.  Otherwise,  we  continue  as  before. 

Now  we  show  that  we  do  not  insert  too 
much  token  weight.  Consider  that  at  any 
level  where  tokens  are  inserted  there  are  less 
than  t--  -t-°'(cns  nodes  in  the  level, 

and  we  need  to  use  at  most  dmax  tokens  to 

push  out  any  node.  Thus  we  use  less  than 
the  number  of  tokens  inserted  tokens  tQ  push 

out  the  nodes.  In  addition,  each  node  that  is 
pushed  out  can  push  out  at  most  dmax  tokens 
with  it  so  at  least  half  the  tokens  that  start  in 
the  level  stay  in  this  level.  Hence,  we  can  add 
only  as  many  tokens  as  we  already  have  which 
means  that  the  number  of  tokens  at  most  dou- 


bles. 

Thus  the  sum  of  the  path  lengths  can  at  most 
double  since  (as  in  Lemma  1)  the  total  distance 
of  path  segments  formed  in  Stage  3  is  0(  + 

s  53,^0  ft,  and  a  ft,  is  upper  bounded  by 
the  number  of  tokens. 

To  see  that  the  edge  depth  of  the  region 
is  0(^*p).  we  restrict  our  attention  to  those 
levels  for  which  the  number  of  nodes  encoun¬ 
tered  exceeds  t r —  of  the  number  of  tokens 
(i.e.,  edges)  encountered.  All  nodes  are  on 
these  levels  by  our  construction.  Since  the 
total  number  of  nodes  reached  decreases  by 
a  factor  of  (1  -  •)  every  time  such  a 

level  is  encountered  during  this  stage,  we  can 
conclude  that  the  edge  depth  of  the  part  of 
the  region  constructed  during  Stage  3  is  at 
most  0 ( n  ) .  Finally  for  bounded  degree 
graphs  this  is  0( =  Of^). 

This  concludes  the  proof  of  Theorem  2. 

4  Remarks 

We  are  in  the  process  of  investigating  possible 
extensions  of  this  work.  Potential  questions  of 
interest  include: 

1.  Can  the  routing  result  of  Section  1.5  be  ex¬ 
panded  to  include  unbounded  degree  net¬ 
works  by  adding  a  diameter  constraint? 

2.  Can  our  methods  be  extended  to  find  ap¬ 
proximate  node  separators? 

3.  Can  the  results  of  Section  1.4  be  extended 

to  approximate  ^  \  separators? 

4.  Can  the  technique  of  Section  1.5  be  ex¬ 
tended  to  improve  the  Peleg-Upfal  disjoint 
paths  result  for  expanders  [15]? 

5.  Are  there  other  combinatorial  problems 
that  can  be  approximated  using  these 
techniques? 

We  suspect  that  many  of  these  questions  can 
be  answered  affirmatively,  but  have  not  yet 
worked  through  the  details. 
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